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I briefly review the problem of representability — a single coarse-grained effective pair potential 
cannot simultaneously represent all the properties of an underlying more complex system - as well 
as a few other subtleties that can arise in interpreting coarse-grained potentials. 
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I. INTRODUCTION 

An important source of progress in computer simula- 
tion comes from better coarse-grained models of the un- 
derlying materials, that is descriptions that are simpler 
and more tractable, but nevertheless retain the funda- 
mental underlying physics that one is interested in in- 
vestigating [lJ-Q ■ To most of us it is intuitively obvious 
that these simplified descriptions throw some information 
away and that compromises are made. After all, there 
is no such thing as a free lunch. But how much are we 
paying, and what can we get away with? 



A. coarse-graining: structure and energy routes 

To illustrate these questions, consider the example of 
a one-component reference fluid interacting with a three- 
body Hamiltonian of the form: 



H = K 



i<3 



i<j<k 



(3) 



(Tij,rj k ,rki), (1) 



where r\ denotes the position of particle i and = r^— Tj 



and 



r, ; - r, 



K is the kinetic energy opera- 



tor, w^ 2 \r) is an isotropic pairwise additive potential, 
and w/ 3 )(rjj, r^, r&j) is a triplet or three-body potential. 
Three-body potentials are expensive and cumbersome to 
simulate, and so one might want to coarse-grain them to 
a simpler isotropic representation. 

There are several ways you could do this. One popular 
method is to fit a pair-potential such that it reproduces 
a structural quantity like the radial distribution function 
g(r) generated by the Hamiltonian of Eq. [TJ Hender- 
son [J] first showed that "the pair potential v(r) which 
gives rise to a radial distribution function g(r) is unique 
up to a constant." . A more rigourous mathematical dis- 
cussion is provided by Chayes, Chayes and Lieb |5|. An 
extended proof for orientational correlations can be found 
in a book by Gray and Gubbins @ and we did this for 
multi-site potentials in ref. [7J. Therefore, for a given 
state-point, the g(r) generated by a Hamiltonian like that 
of Eq. ([T]) can be reproduced by a unique effective pair 
potential v g (r). (The existence of v g is more subtle, but 
holds under fairly general conditions, see @). I'll call ap- 
proaches that derive v g (r) the structural route to deriving 



an effective potential. The difference with the bare-pair 
potential (r) can be written as: 



Sv g (r) 



,(2) 



(r) -v g (r). 



(2) 



Another popular way of deriving an effective poten- 
tial is by fitting to a thermodynamic observable like the 
internal energy. One first calculates the full internal en- 
ergy U(N,V,T) for the system governed by the Hamil- 
tonian (TT|), and then (in this case) derives a potential 
v u^ ( r ) that reproduces the same internal energy by the 
simpler two-body formula: 



U(N,V,T) = l -p 2 J dv x dv 2 g{r 12 )vf{r l2 ). 



(3) 



I'll call this the energy route. The difference with the 
bare-pair potential w^ 2 \r) can be written as: 



5vu(r) 



= ,„(3) 



(r) - vu(r). 



(4) 



It is not hard to show that both v g (r) and vjj (r) depend 
on the state point at which they are derived. Thus if 
they are used at a different state-point, one would expect 
transferability problems, i.e. you'd need to re-derive 
them for your new state-point. 



B. representability problems 

What is perhaps more worrying is that one can also 
show that at a given state-point, v g (r) and vu{r) can- 
not be the same. As first demonstrated almost 40 years 
ago [8[, to lowest order in p and u/ 3 \ the ratio between 
the two corrections is: 



Svu(r) 
Sv g (r) 



O 



(^ (3) ) 2 ;p 2 



(5) 



In other words even in the low density weak potential 
limit, the corrections due to the three-body forces that 
you are coarse-graining out differ by a factor of three! 
What happens for stronger interactions or higher densi- 
ties needs to be investigated. Since v g (r) is unique, it is 
therefore impossible to represent all the properties of a 
system governed by the Hamiltonian of Eq. ((T|) by a sin- 
gle pair potential. There is no free lunch. I believe such 
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representability problems are widespread in coarse- 
grained descriptions of soft- matter systems [9(. 

They may also be important important in other coarse- 
grained simulations. For example, we recently used the 
structural route to derive radially symmetric potentials 
for water from a more sophisticated underlying model 0] ■ 
We explicitly constructed v g (r) and vu(r) and they look 
very different, as anticipated by Eq. [51 At some of the 
state-points we studied, using v g (r) to calculate the virial 
pressure resulted in a dimensionless compressibility fac- 
tor Z = fiP/ p that was almost two orders of magnitude 
larger than that of the original multi-site water model 
used to parameterise v g (r). Similar deviations from the 
underlying water model were also seen in ref. (ip| . 

Admittedly, it may not be surprising that an isotropic 
potential should perform so poorly when the underly- 
ing fluid has complex oricntational correlations. In fact 
treating water as a pair potential is probably not such a 
good idea 

One way of improving on the thermodynamic perfor- 
mance of the structurally derived potential v g (r) is to 
use constraints to simultaneously fit to properties like 
the virial pressure or the internal energy. Due to the 
uniqueness of v g (r), this potential would no longer cor- 
rectly reproduce the pair correlations, but since the ther- 
modynamic properties are scalar quantities, it might not 
come at too large a cost for the structure. A nice exam- 
ple of applying this to water can be found in this paper 
by the Mainz group who showed that the structural 
properties were affected (and by extension the isothermal 
compressibility which can be written as an integral over 
g(r)), but not too badly. 



C. more general problems with effective potentials 

In ref. @ some more general issues with effective po- 
tentials are reviewed. For example, most of the time 
our effective simplified potentials are really state depen- 
dent, that is if you were to derive them at different state 
points they would be different. I argue that you shouldn't 
naively treat them as if they are real two-body potentials 
of a Hamiltonian system. 

Here is an example of how things can go horribly 
wrong. Consider a homogeneous fluid in a volume V, 
with N particles interacting with a spherically symmet- 
ric pair potential v(r; p), that depends on the state point 
(that could be density, temperature, etc.). For simplic- 
ity, here we consider only the density p = N/V. The 
standard way to derive the virial equation is directly 
through the canonical partition function 



Q(N,V,T) 



A 



-3N 



J dr N exp |-^^w(r li; p)| 



volume derivative in 



(3P = 



d log Q(N,V,T) 
dV 



(7) 



N,T 



should also act on v(r;p), resulting in a virial equation 
with an extra dv(r; p)/dp term: 



7 P 

vir 



2 

JQ 



dv(r;p) 
dr 



(8) 

P dp J g(r) ' 



( 6 ) 

where A is the usual thermal de Broglie wavelength. The 



something first pointed out in 1969 by Ascarelli and 
Harrison [ll| in the context of density dependent pair po- 
tentials used for modelling liquid metals. 

Now on the surface this all looks very kosher - I 
took my potential, plugged it into my partition function, 
turned the statistical mechanical crank, and out popped 
Eq. \E\ with a correction to the pressure due to the state- 
dependence of the effective potential. However, when we 
tried this with the v g (r) derived for the water model Q, 
the correction actually made the agreement far worse. 
In ref. Q I give other examples were this correction has 
the wrong sign (and magnitude). Although the analysis 
leading to Eq. [5] has the veneer of statistical mechanical 
respectability, it is in fact a disreputable result. 

Taking a step back, it is not at all surprising that this 
derivation cannot be right The correction term in 
Eq. [5] only takes into account a local dependence on den- 
sity. If you took the same potential v(r; p), and used it in 
a grand-canonical ensemble, then you would need infor- 
mation from all densities p, not just the local derivative. 
So the two ensembles would not be equivalent, a good 
sign that there are more problems afoot. Warnings about 
such difficulties abound, for example, about 40 years ago 
Barker et al. [l2j wrote: 

We record our opinion that the use of 
density- dependent effective pair potentials can 
be misleading unless it is recognised that these 
are mathematical constructs to be used in 
specified equations rather than physical quan- 
tities 

and similar reservations were sounded by other distin- 
guished investigators @, HH, EH • 

The point of this exercise is not to say that effective po- 
tentials are useless. I've happily used them myself (la ]. 
Sometimes a fine solution to a dubious Hamiltonian is 
better than a dubious solution to a very fine Hamilto- 
nian. It is just that one must be careful when using 
coarse-graining procedures not to automatically treat the 
coarse-grained potentials as if they are real bona-fide po- 
tentials of the type one uses in a Hamiltonian. Also, in 
practice it is also good to remember that these potentials 
are always compromises and that fitting too strongly to 
one quantity may generate larger errors in other quanti- 
ties that one also wants to measure 0, ■ 
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II. QUESTIONS FOR ALL 

1. Do others have good examples of representability 
problems generated by their potentials? I'd be in- 
terested to know how important they are in other 
systems. 

2. A related question is: How well can you predict 
the magnitude of a representability problem? In 
the case of water it is not surprising that a radially 
symmetric pair-potential would have problems, but 



could you have predicted that a potential that gen- 
erates g(r) so well does so poorly on the internal 
energy or the virial pressure? 

3. Ignacio Pagonabarraga has an interesting interpre- 
tation of density dependent potentials in ref. [16j . 
Do you think this approach could be generalised to 
the kind of coarse-graining procedures we are ap- 
plying in this conference? 
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